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The  aim  of  the  present  work  is  the  calculation  of  the  energy  of  ionic 
xenon-chlorine  systems  which  can  be  formed  in  solid  Xe  by  irradiation.  The 
energy  levels  of  these  ionic  systems  differ  from  those  in  the  gas  phase  due 
to  polarization  aa^dispersion  interactions  with  solid  Xe  atoms.  It  is 
shown  that  the  Ke*Cl  y)>olecule  is  responsible  for  experimentally-observed 
emission.  The  activation  energy  of  thefxe^Cl  formation  is  found  to  form  a 
broad  band. 
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I.  Introduction 

Excited  ionic  rare  gas-halogen  molecules  are  of  great  interest, 
particularly  because  of  their  lasing  properties.  These  molecules  have  been 
well  investigated  in  the  gas  phase,  both  experimentally  and 
theoretically.  *  ^  Recent  experiments  performed  by  Fajardo  and  Apkarian  ^ 
have  demonstrated  that  ionic  xenon-chlorine  molecules  can  also  be  formed  in 
a  Xe  solid  matrix.  The  spectral  transitions  in  these  molecules  in  the  solid 
phase  are  not  the  same  as  in  the  gas  phase  due  to  interactions  between  the 
ionic  molecules  and  the  surrounding  solid  atoms.  Such  interactions  are 
considered  in  the  present  paper  in  order  to  study  the  electronic  transitions 
in  the  ionic  xenon-chlorine  molecules  inside  the  Xe  solid  matrix. 

The  paper  begins  with  the  experimental  background  in  Section  II,  both 
gas  phase  and  solid  phase,  and  then  proceeds  with  the  theoretical  analysis 
of  the  molecule-matrix  interaction  in  Section  III.  The  numerical  analysis 
of  the  Xe^Cl  molecule  and  the  activation  complex  for  ionic  state  formation 
is  presented  in  Sections  IV  and  V,  and  Section  VI  is  the  Conclusion. 

II.  Experimental  Background 

A.  Xe+Cl  and  Xe^Cl  Molecules  in  the  Gas  Phase 

The  formation  of  excited  molecules  (XeCl)  and  (Xe^Cl)  by  laser 

1-3 

irradiation  has  been  observed  in  Xe/C^  gas  mixtures.  These  molecules 

are  ionic  in  nature  and  hence  are  denoted  by  Xe+Cl  and  Xe^Cl  . ^  ^  The 

Xe+Cl  molecule  is  formed  as  a  result  of  the  absorption  of  two  separate 

ultraviolet  photons.  The  first  of  these  photons  hits  a  Cl^  molecule  and 

produces  by  photodissociaton  the  Cl  atoms.  The  second  photon  at  hv  *  4.02 

eV  provided  specifically  by  an  XeCl  excimer  laser  is  absorbed  by  a  colliding 

Xe  +  Cl  pair  transferring  an  electron  from  Xe  to  Cl  and  forming  the  excited 
+  -  2 

Xe  Cl  molecule.  The  electronic  structure  of  the  XeCl  system  is  shown  in 


Fig.  1.  The  two  lowest  states,  X  and  A,  are  formed  by  neutral  Xe  and 
2 

C 1  (  P)  atoms.  The  ground  state  is  weakly  bound  by  0.035  eV  and  becomes 

°  8 

repulsive  at  relatively  short  interatomic  distances,  R  <  3.1  A.  The  two 
excited  ionic  states  are  similar:  both  have  a  deep  minimum  and  form  excited 
ionic  molecules  due  to  the  attractive  electrostatic  forces.  The  photon 
excites  the  system  to  the  B  state,  but  due  to  the  crossing  of  potential 
curves,  the  molecule  can  go  to  the  excited  C  state  from  which  a  radiative 
transition  to  the  repulsive  A  state  becomes  possible  with  the  emission  of  a 
photon  of  3.54  eV  (X  =  350  nm). 

The  emission  of  the  350  nm  photons  is  partly  quenched  by  the  formation 

of  the  triatomic  ionic  molecules  Xe^Cl  due  to  the  three-body  collisions  of 

the  Xe+Cl  molecule  with  two  Xe  atoms. ^  The  Xe^Cl  molecule  is  stabilized 

relative  to  the  Xe+Cl  molecule  by  a  strong  Xe+-Xe  valence  attraction. 

2 

Due  to  the  lower  energy  of  the  excited  4  r  state,  radiative  emission  from 

^  ^  “ 

Xe^Cl  is  red-shifted  compared  to  the  Xe  Cl  emission.  The  Xe2Cl  emission 

is  observed  as  a  wide  band  with  its  center  at  A  =  480  nm  (2.58  eV), 

[Xe^Cl* ]^2p  ->  Xe  +  Xe  +  Cl  +  hv(-2.58  eV)  .  (1) 

Energy  levels  of  the  Xe^Cl  system  have  been  calculated  by  the 
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diatomics- in-molecules  method,  ’  and  the  potential  curves  are  shown  in 

2 

Fig.  2.  The  equilibrium  point  of  the  lower  ionic  state  (4  r)  was  obtained 

o  o 

for  the  interatomic  distances  *  3.32  A  and  ^e_^e  =  3.29  A,  where  the 

three  atoms  form  a  triangle  with  almost  equal  sides.  The  calculated  Xe-Cl 


interatomic  distance  is  larger  than  the  interatomic  distance  in  the  diatomic 

+  -  0 
Xe  Cl  molecule,  where  the  latter  has  been  experimentally  measured  as  2.94  A 


°  4  + 

and  calculated  as  3.22  A.  However,  the  Xe-Xe  distance  in  the  Xe2Cl 

molecule  is  close  to  the  interatomic  distance  of  the  Xe*  molecule, 


3.27  A  . 


The  Xe0Cl  molecules  demonstrate  lasing  at  X  *  518  nm,  which  is 
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na. 
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B.  Formation  of  Ionic  Xenon-Chlorine  Emission  Centers  in  Xe  Solid 

The  formation  and  emission  of  the  ionic  xenon-chlorine  centers  has  been 
observed  in  Xe  solid  by  Fajardo  and  Apkarian.^  The  experiments  were 
performed  with  Xe  solid  containing  molecules  as  an  impurity  in  the 
proportion  1/100-1/5000.  After  the  solid  was  irradiated  by  308  nm  laser 
light,  the  emission  was  detected  at  A  *  570  nm.  This  emission  can  be 
interpreted,  as  in  the  gas  phase,  by  the  photodissociation  of  followed 
by  the  formation  of  the  excited  ionic  Xe-Cl  system  due  to  the  absorption  of 
the  next  photon.  According  to  the  experiment,^  only  one  excited  ionic 
center  (most  probably  Xe^Cl  )  is  formed  since  only  one  emission  band  (A  = 

570  nm)  is  observed  in  the  solid,  contrary  to  the  gas  mixture  where  two 
ionic  molecules  (Xe+Cl  and  Xe^Cl  )  can  be  formed  with  two  different 
emission  wavelengths.  The  Xe+Cl  molecules,  if  they  are  formed  in  the 
solid,  are  expected  to  relax  quickly  to  the  more  stable  Xe^Cl  molecule. 

The  experiments  have  also  been  performed  with  Xe  solid  containing  HC1 
molecules,  and  these  demonstrate  the  same  spectroscopic  properties  as  in  the 
case  of  CI2  +  Xe  solid  system.  The  conclusion  can  thus  be  made  that  after 
the  photodissociation  of  either  or  HC1,  the  atoms  of  these  molecules  are 
completely  separated  in  the  solid  without  any  influence  on  one  another. 
Simultaneously  with  the  optical  emission,  the  +  Xe  and  HC1  +  Xe  solid 
systems  demonstrate  electrical  conductivity,  implying  that  the  formation  of 
the  excited  ionic  centers  is  accompanied  by  the  excitation  of  current 
carriers.^ 

We  shall  consider  in  the  present  work  below  the  following  problems: 

(1)  the  calculation  of  the  wavelength  shift  in  the  Xe^Cl  emission  in  the 


solid  compared  to  the  gas-phase  emission;  and  (2)  the  consideration  of  the 
activation  complex  and  the  calculation  of  the  excitation  energy  for  the 
Cl  +  Xe  solid  system  transition  to  the  ionic  state.  We  shall  not  take  into 
account  the  vibrational  motion  of  the  solid  atoms. 


III.  Molecule-Rare  Gas  Solid  Matrix  Interaction 

By  a  molecule  in  a  rare  gas  solid  matrix,  we  understand  a  system  of 
atoms  (ions)  with  much  stronger  inner  interactions  than  the  interactions 
between  the  molecule  and  the  matrix  atoms.  This  definition  is  obvious  for 
the  case  of  neutral  atoms  or  strongly  bonded  molecules,  like  or  HC1 
doped  in  rare  gas  solid.  However,  for  the  case  where  the  molecule  includes 
rare  gas  atoms  of  the  host  lattice  (Xe+Cl  molecule,  for  example),  the 
definition  of  the  molecule  may  be  ambiguous. 

Considering  molecules  inside  a  rare  gas  matrix  we  shall  make  three 
important  assumptions:  (i)  the  electronic  bands  of  the  rare  gas  solid  are  so 
narrow  that  the  solid  structure  does  not  change  the  interatomic  potential 
curves  from  those  in  the  gas  phase; ^  (ii)  the  valence  (exchange) 
interactions  between  the  atoms  of  the  molecule  and  the  matrix  atoms  can  be 
neglected;  and  (iii)  the  coupling  between  the  neutral  and  ionic  states  of 
the  molecule  can  be  neglected,  and  therefore  the  neutral  and  ionic  molecules 
(Xe2Cl  and  Xe^Cl  ,  for  example)  can  be  considered  independently.  The  first 
of  these  assumptions  is  well  defined,  so  we  shall  discuss  the  two  other 
assumptions  only. 

The  valence  interaction  between  a  molecule  and  matrix  can  be  important 
in  the  case  of  ionic  molecules  R*  and  R*X  (X  and  R  are  halogen  and  rare  gas 
atoms,  respectively)  where  the  positive  charge  is  partly  delocalized  between 
the  molecular  rare  gas  atoms  R+  and  the  matrix  atoms.  This  valence 
interaction  has  been  discussed  in  the  literature  for  the  case  of  R^  molecule 
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in  an  R  matrix,  and  in  Ref.  13  the  neglect  of  the  valence  interaction 

was  proposed.  The  error  resulting  from  the  neglect  of  the  molecule-matrix 

valence  interaction  for  Xe+Cl  molecule  in  Xe  matrix  will  be  estimated  in 

n 

Section  V. 

The  coupling  between  the  neutral  and  ionic  states  of  the  trapped 

molecule  is  weak  due  to  the  high  energy  of  the  rare  gas  atoms  excitation  and 

small  overlap  of  orbitals.  The  weakness  of  the  coupling  in  the  Xe-halogen 

atom  interaction  is  well  demonstrated  by  the  results  cf  an  ab  initio 
4 

calculation.  This  calculation  examines  the  charge  transfer  between  the  Xe 
and  halogen  atoms  which  can  be  considered  as  a  measure  of  the  coupling.  The 
charge  transfer  decreases  with  the  interatomic  distance,  and  for  the  Xe-Cl 
system,  for  example,  it  becomes  smaller  than  10 7.  when  the  Xe-Cl  distance 

O 

exceeds  3  A.  In  the  Xe  matrix  the  distance  between  trapped  atoms  and  Xe 

O 

atoms  is  usually  larger  than  3  A,  and  consequently  the  coupling  is  small. 

The  weakness  of  the  coupling  is  supported  also  by  an  experimental  study  of 

the  spectrum  of  homonuclear  halogens  trapped  in  a  rare  gas  matrix.  The 
(s) 

interaction  U  between  these  molecules  and  the  solid  matrix  is  mainly  of 
AM 

van  der  Waals  origin,  and  therefore  it  is  small  and  does  not  differ 
significantly  for  different  electronic  states.  It  follows  that  in  the 
absence  of  coupling  the  molecular  electronic  states  are  expected  to  be 
practically  the  same  as  in  the  gas  phase,  in  accord  with  the  experiment.^ 
The  above  assumptions  greatly  simplify  the  consideration  of  the  system, 
since  the  total  energy  of  the  molecule  (A)  and  the  matrix  (M)  can  be 
expressed  as  a  sum  of  the  molecular  energy  the  solid  matrix  energy 

(E^S^)  and  the  energy  of  the  molecule-matrix  interaction 


i?(s)  4.  ij(s) 


C 


rv 


The  molecular  energy  E^®  is  the  same  as  in  the  gas  phase  for  the 

corresponding  geometry,  which  can  be  different  than  the  equilibrium  gas 

(s) 

phase  geometry.  The  matrix  energy  is  the  energy  of  interaction  between 

matrix  atoms  excluding  those  of  them  which  belong  to  the  molecule.  In  Eq. 

(2)  and  below  the  values  associated  with  the  gas  phase  and  solid  state  are 

distinguished  by  the  upper  indexes  (g)  and  (s),  respectively. 

For  the  case  of  a  neutral  molecule  which  does  not  have  a  dipole  moment, 
(s) 

the  term  in  Eq.  (2)  is  equal  to  the  sum  of  diatomic  interaction 

AM 

potentials  between  the  molecule  and  the  matrix  atoms.  For  the  case  of  a 
molecule  which  generates  an  electric  field,  in  particular  an  ionic  molecule, 
the  polarization  of  the  matrix  atoms  has  to  be  taken  into  account.  The 
polarization  energy  can  not  be  presented  as  a  sum  of  diatomic  interactions, 
so  it  is  calculated  separately  whereas  the  polarization  terms  are  excluded 
from  the  diatomic  potentials: 


i  J 


Here  is  the  polarization  energy,  i  and  j  are  atomic  indices  of  the 
molecule  and  the  matrix,  respectively,  and  are  the  interatomic 
interactions  without  polarization  terms, 


°ij  *  Uu  •  '-V*  u)> 


u> 


where  U. .  is  the  diatomic  potential,  R, .  is  the  interatomic  distance,  and 
ij  r  ij 

-A  /R^ .  is  the  polarization  term  which  is  valid  for  ion-atom  and  ion-ion 
P  ij 

interactions  only.  When  both  atoms  are  neutral,  0,.  coincides  with  U.  . . 

ij  ij 
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The  energy  of  the  solid  matrix  polarization  can  be  calculated  by 
forning  around  the  ionic  molecule  an  Onsager  cavity  and  considering  the 
solid  matrix  outside  the  cavity  as  a  continuum,^ 

e  ■  e  +  e  ,  (5) 

p  p,c  p,o 

where  ^  is  the  polarization  energy  of  the  atoms  inside  the  cavity  and 

e  is  the  polarization  energy  of  the  continuum  outside  the  cavity.  The 
p,o 

dipole  moments  of  the  polarized  atoms  are  assumed  to  be  formed  by  charged 
atoms.  The  polarization  energy  of  the  atoms  inside  the  cavity  includes  also 
the  dipole-dipole  interaction 


=-A  Jz2  +  2aA  y  y  Z.  Z.  , 

>,c  p  L  L  J*m  P  L  L  L  L  ij.mn  i.m  J.n 

j=l  m  i  j ( >i )  m  n 


where  m  *  1,2,3  denote  the  x,y,z  projections  of  the  2^  vectors  which  are 
proportional  to  the  field.  When  the  field  is  formed  by  a  diatomic  ionic 
molecule,  then  J  *  2  and 


5.-S. 

t. - Vs  +  -J-r  ' 

1  R?, 


When  the  field  is  formed  by  a  triatomic  ionic  molecule,  like  Xe2Cl  ,  then 
J-3  and 

2  =  .  _j _ 0  +  _j _ 1  +  _j _ 2  .  w  ) 

j  R» .  2R? .  2R^ . 

Oj  lj  2j 


The  matrix  of  the  dipole-dipole  interaction  is 


v  ■»■»  ir»  f  -r-  r»  r»  ^  '.■vy*-'t^  'A  a  '.v.r.r.v.v.rT.v.r  v-.v.^-w  v  v  vwv\vv  -.- «.-  v  ■ 


Y,  .  =  3($.-  J.)  (t.  -  1)  /R5.  -  6  /R3., 

'ij,mn  j  i  m  j  i  n  ij  mn  ij 


(8) 


£ 


s\- 


where  6  is  the  Kronecker  delta.  In  the  expressions  ( 6 ) - ( 8 )  i  and  j  are 


mn 


the  indices  of  the  matrix  atoms,  is  the  position  vector  of  the  negative 
ion,  and  and  ^  are  t^ie  position  vectors  of  the  positive  rare  gas  ions. 
The  parameter  is  proportional  to  the  rare  gas  atom  polarizability  a. 


A  =  ae  /  2 . 
P 


(9) 


The  contribution  of  the  dipole-dipole  interaction  to  the  polarization  energy 

18 

(6)  depends  on  the  atoms  localization  and  the  dielectric  constant.  In  the 
present  calculation  the  contribution  of  the  dipole-dipole  interaction  stands 
for  25-30  7.  of  the  total  polarization  energy.  The  polarization  of  the 


continuum  outside  the  cavity  is  considered  in  the  dipole  approximation. 
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P.o 


i  2  2 

_  e  ~  1  e  P 
2c  +  !  r3 


(10) 


where  e  is  the  dielectric  constant  of  the  solid,  p  is  the  distance  between 
opposite  charges  within  the  ionic  molecule,  and  is  the  radius  of  the 
cavity. 

The  atom-atom  potentials  are  expressed  in  this  paper  either  by  the 
Lennard- Jones  expression 

U(R)  «  cQ  [(Rq/R)12  -  (Rq/R)6]  (11) 

or  by  representing  the  repulsive  term  by  an  exponential  function, 

U(R)  -  A  exp(-aR)  -  A./R6,  (12) 

r  d 


9 


where  A,  is  the  dispersion  coefficient.  The  ion-atom  potential  includes 


also  the  polarization  term 


U(R)  =  A  exp(-aR)  -  A./R6  -  A  /R4.  (12') 

r  a  p 

The  polarization  coefficient  is  determined  by  Eq.  (9).  The  energies  of  the 
electronic  states  are  denoted  below  by  E,  the  interatomic  potentials  by  U, 
the  energies  of  system  formation  by  V,  and  the  transition  energies  by  T. 

TV.  Xetci  Molecule  in  a  Xe  Solid  Matrix 

The  crystal  structure  of  the  Xe  solid  is  face-centered  cubic  (fee)  with  four 

atoms  per  unit  cell.  Each  atom  in  the  lattice  is  surrounded  by  12  nearest 

neighbors  with  the  distance  2.  =  a//2  and  6  next  neighbors  with  the  distance 

19 

a,  where  a  is  the  lattice  constant.  The  Xe  lattice  constant  is 

O 

a  =  6.197  A  for  T  ~  58  K,  and  correspondingly  the  minimal  interatomic 

°  20 

distance  is  l  =  4.38  A.  The  minimal  distance  t  extrapolated  to  0  K  is 

o  2\ 

only  0.7%  smaller  (4.35  A).  All  further  calculations  are  performed  with 
i  =  4.38  A. 

°  22 

The  radius  of  a  negative  ion  Cl  is  1.82  A,  -83%  of  the  Xe  atom 

o  _ 

radius  (2.19  A),  so  it  is  possible  to  assume  that  Cl  replaces  one  Xe  atom 
without  deforming  much  the  surrounding  lattice.  Because  of  the  small 
difference  between  the  radii  of  Cl  and  Xe,  the  Cl  ion  can  be  shifted 
inside  the  lattice  cage  by  small  distances  only.  We  thus  assume  that  Cl  is 
located  at  a  lattice  point,  whereas  the  relaxation  of  the  Xe  atoms  is 
admitted.  The  formation  of  the  Xe^Cl  molecules  has  to  be  considered  now  as 
a  shift  of  two  neighboring  Xe  atoms  with  a  common  charge  of  +1  towards  the 
Cl  ion. 

The  interaction  energy  between  the  atoms  inside  the  Xe^  Cl  molecule  is 
much  larger  than  the  interaction  energy  between  the  Xe+Cl  molecule  and  the 
Xe  matrix,  so  it  is  reasonable  to  assume  that  the  geometry  of  the  Xe^Cl 


molecule  in  the  solid  matrix  is  the  same  as  in  the  gas  phase.  Calculational 

results  support  this  assumption  since  the  interatomic  distances  of  the 

Xe^Cl  molecule  in  solid  are  found  to  be  larger  than  those  in  the  gas  phase 
°  +  - 

by  0.1  A  only.  As  the  Xe^Cl  molecule  has  the  form  of  a  triangle  with 
almost  equal  sides  where  the  angles  are  close  to  60°,  the  two  Xe  atoms  which 
form  the  Xe^Cl  molecule  are  shifted  in  the  lattice  along  crystallographic 
lines.  The  coordinates  of  some  of  the  atoms  are  presented  in  Table  1.  The 
Cl  ion  is  located  at  the  origin  of  the  coordinate  system  *  0). 


A •  Energy  of  the  Xe^Cl  -*•  Xe^Cl  Transition 

The  energy  of  the  transition  (1)  in  the  Xe  solid  is 

(s)  (s) 

^e^Cl*  =  ^e^Cl 

(s) 


E,(S) 

T(e2Cl  ' 


(13) 


(s) 


Here  E^e+C^-  is  the  Xe2Cl  energy  in  solid,  and  E^g  ^  is  the  energy  of  two 


neutral  Xe  atoms  and  one  Cl  atom  in  'be  solid  with  the  interatomic  distances 


equal,  according  to  the  Franck-Condon  approximation,  to  the  corresponding 


distances  in  the  Xe^Cl  molecule.  The  interatomic  interactions  0„  in  Eq. 


(3)  can  be  restricted  to  the  nearest-neighbor  pairs.  For  the  Xe^Cl  case 
they  include  the  interactions  of  the  Cl  ion  with  10  Xe  atoms  located  around 


the  ion  and  the  interaction  of  each  of  two  Xe 
matrix  atoms, 


1/2+ 


atoms  with  10  nearest 


(s)  (*>  (s) 

^Cl*  =  ■x.Jcr  +  "m  +  10  °XeCl" 


Xe2M 


+  UXe+M  +  Gp  ' 

(14) 

'l,16)  ’ 

(15) 

1/2+ 


where  0..  stands  for  the  interaction  between  the  Xe  atom  (i  =  1)  and 
ij 


different  matrix  atoms  with  different  interatomic  distances  (see  Table  1), 


12 


The  Xe^Cl  energy  is 


p(s) 

Xe^Cl 


I?(  8  )  a. 

Ae2Cl 


E*s)  +  10U 


XeCl  +  2UXeH  ’ 


(16) 
),  (17) 


UXeM  "  3UXeXe(Rl,3)  +  UXeXe(Rl , 13)  +  4UXeXe(Rl , 14J  +  2UXeXe(Rl,16 
(e) 

where  E£*  ^  is  the  energy  of  the  Xe2Cl  system  in  the  configuration  of  the 

Xe^Cl  molecule,  and  the  R. .  are  the  distances  between  the  i  *  1  Xe  atom  and 
2  lj 

nearest-neighbor  matrix  atoms  (see  Table  1).  In  accord  with  the  Franck- 

(s) 

Condon  approximation,  the  matrix  energy  E^  is  the  same  in  Eqs.(14)  and 
(16). 

In  order  to  calculate  the  energy  of  transition  (13)  one  needs  to  know 
the  polarization  energy  and  the  following  diatomic  potential  curves: 
l'x  ci»  Oxed'*  UXeXe*  and  ^Xe^2+Xe‘  The  Polarization  energy  (5)  will  be 

o 

calculated  using  the  radius  of  the  Onsager  cavity  r^  =  7.83  A.  The  cavity 

with  this  radius  includes  centers  of  all  23  (i  M  3, 4,..., 25)  Xe  atoms  which 

are  the  nearest  neighbors  of  the  atoms  of  the  Xe2Cl  molecule,  and  only  one 

of  these  is  close  to  the  cavity  border.  Approximately  8  atoms  located 

wholly  or  partly  inside  the  cavity  are  not  included  in  the  sum  (6),  since 

these  atoms  are  not  the  nearest  neighbors  of  the  Xe^Cl  molecule,  and  thus 

their  contribution  to  the  polarization  energy  is  small  (less  than  10%).  The 

20 

dielectric  constant  of  Xe  solid  is  e  *  2.1,  and  the  distance  between 

+  -  8 

opposite  charges  within  the  Xe2Cl  molecule  is  p  ■  2.86  A.  The  polarization 

coefficient  (9),  which  is  proportional  to  the  Xe  atom  polarizability  a  * 

°3  23  °4 

4.05  A  ,  is  equal  to  A^  *  29.14  eV  A  .  Substituting  A^,  e,  p  and  Tq  into 

the  expressions  (6)-(8),  one  obtains  e  ■  -0.34  eV,  t  *  -0.05  eV  and  e  * 

p,c  p,o  p 

-0.39  eV. 

The  Xe-Cl  potential  is  expressed  by  Eq.  (12)  which  fits  the  ground 

O 

state  potential  of  Ref.  8  for  R  >  3  A.  The  coefficients  of  the  Xe-Cl 


1 


iVtt.'iUo. 


°-l  °6 

potential  are  found  to  be  a  *  2.3  A  »  Af  *  2A8  eV  and  A^  3  206  eV  A  .  The 

24 

experimental  potential  curve  of  the  Xe-Xe  interaction  is  well  fitted  by 

the  Lennard- Jones  potential  (11)  with  the  coefficients  3  0.096  eV  and  Rq 
*  3.98  A. 

The  halogen  ion-rare  gas  atom  interaction  potential  was  suggested  to  be 

expressed  as  a  sum  of  a  repulsive  exponential  term  and  an  electrostatic 

(polarization)  term.  ’  However,  these  two  terms  do  not  fit  well  the 

known  halogen  ion-rare  gas  atom  potentials,  so  that  the  expression  (12') 

with  a  dispersion  term  was  used  to  fit  the  recent  experimental  Xe-Cl 

potential.  The  coefficients  of  this  potential  were  found  to  be  a  -  3.85 
0  1  0£ 

A  ,  A  *  31670  eV  and  A  =  60  eV  A  .  The  polarization  coefficient  is  A  = 
r  d  p 

29.14  eV  A  .  The  discrepancy  between  the  analytical  potential  (12')  and  the 

O 

experimental  values  does  not  exceed  0.005  eV  for  R  >  3.1  A.  The  potentials 

28 

fcr  some  other  halogen  ion-rare  gas  atom  pairs  can  be  also  fitted  by  the 

potential  (12’).  The  Xe-Cl  potential  0  which  appears  in  the  Eq.  (14)  is 

determined  by  the  expression  (12),  i.e.,  without  the  polarization  term  - 

A  /R4. 

P 

1/2+ 

Unfortunately,  the  Xe  -Xe  U  potential  cannot  be  determined  in  the 
same  way  as  in  the  Xe-Cl  case.  The  Xe+-Xe  potential  is  known, ^  but  the 
nature  of  this  potential  is  complicated  and  includes  valence,  electrostatic 
polarization  and  dispersion  effects  which  cannot  be  separated  from  one 
another.  However,  taking  into  account  that  the  size  of  the  Xe  ion  is  not 
supposedly  much  different  from  that  of  the  Xe  atom,  we  assume  that  the 
repulsive  and  van  der  Uaals  interactions  between  Xe+  and  Xe  (or  Xe^^  and 
Xe)  atoms  is  the  same  as  between  two  neutral  Xe  atoms.  According  to  its 
definition  [Eq.  (4)],  the  0  potential  of  the  Xe+-Xe  interaction  does  not 
include  the  polarization  term.  It  does  not  include  also  the  valence 


/  j 


attraction  due  to  the  neglect  of  charge  delocalization  between  a  molecule 
and  the  matrix.  It  follows  that  the  repulsive  and  van  der  Waals  terms 
remain  only.  Assuming  that  these  terms  are  equal  to  the  corresponding  terms 
of  the  Xe-Xe  interaction,  one  obtains 

°Xe+X.  ’  °Xe1/2+Xe  *  UXeXe'  (18) 

1/2+ 

The  equality  of  the  Xe  -Xe  0  potential  to  the  Xe-Xe  potential 
greatly  simplifies  the  expression  (13)  for  the  energy  of  transition. 
Substituting  in  (13)  the  expressions  (14)  and  (16),  we  obtain  the  transition 
energy  expression,  which  does  not  include  any  Xe-Xe  interactions, 

(s)  (g) 

TXe*Cr  =  TXe*  Cl"  +  ep  +  1O[0XeCr(8,)  '  UXeCl(e,)1  ’  (19) 

(g)  +  . 

where  Tv  +  =  2.56  eV  is  the  transition  energy  in  the  Xe_Cl  gas 

Xe^  Cl  *■ 

2  o 

molecule  and  l  =  4.38  A  is  the  minimal  interatomic  distance  in  the  Xe 
lattice.  It  is  assumed  in  Eq.  (19)  that  all  Xe  matrix  atoms  are  located  at 
the  lattice  points.  This  assumption  is  supported  by  the  results  of 

O 

calculations  which  indicate  a  small  shift  of  -0.1  A  only  for  atoms 

surrounding  the  Xe^Cl  molecule.  Substituting  into  the  expression  (19),  the 

values  of  the  polarization  energy  (e^  ■  -0.39  eV)  and  the  potentials  (O^cCl* 

(O  «  -0.006  eV,  Uv  _,(fc)  =  -0.019  eV),  one  obtains  that  the  transition 
XeCl 

energy  of  Xe^Cl  in  the  solid  is  T^g  +^-  =  2.30  eV  (Fig.  3),  and 
consequently  the  wavelength  is  A  =  540  nm.  This  value  is  close  to  the 


experimental  wavelength  of  the  center  of  the  emission  band  (570  nm).  The 
shift  in  the  transition  energy  within  the  solid  relative  to  the  gas-phase 
value  is  -0.26  eV  as  compared  to  the  experimental  value  of  -0.38  eV. 


We  now  calculate  the  energy  of  the  formation  of  the  excited  Xe2Cl 
state  from  the  initial  Cl  +  Xe  solid  system,  with  a  Cl  atom  in  the  center  of 
a  cage.  The  calculation  shows  that  in  the  initial  system  the  relaxation  of 

O 

12  Xe  atoms  of  the  cage  towards  the  Cl  atoms  is  small, i.e.,  only  0.06  A. 

This  relaxation  is  neglected,  and  consequently  all  Xe  atoms  are  located  at 
lattice  points.  In  the  initial  Cl  +  Xe  solid  system  we  consider  separately 
a  neutral  Xe2Cl  system  with  the  energy 

-  eci  +  2Exe  +  >4S)  +  i2wi)  +  2uwi)  •  <2o) 

The  energy  of  formation  is  equal  to  the  difference  between  the  Xe^Cl  energy 
(14)  and  the  Xe2Cl  energy  (20), 

(s) 

^Xe^Cl"  =  ^Xe^Cl”  "  EXe2Cl(l!') 

(g) 

"  VXe^Cl-  +  1O0XeCr(*>  ‘  l2UXeCl(i)  +  UXe+M  ‘  21UXeXe(t)  +  Ep  •  (21) 

where  Uv  +u  is  determined  by  expression  (15)  and 
Xe2n 

(g)  (g) 

VXe^r  *  ^*01-  -  (EC1  +  2EXe>  <22) 

is  the  energy  of  the  formation  of  the  Xe+Cl  gas  molecule,  equal  to  3.62 

10  0 
eV.  The  potentials  in  the  expression  (20)  for  t  *  4.38  A  are:  * 

-0.006  eV,  Uv  *  -0.019  eV,  Uv  +u  =  -0.24  eV,  Uv  v  *  -0.024  eV. 

XeCl  Xe2M  XeXe 

(s) 

Inserting  these  values  and  in  Eq.  (21),  one  obtains  =  3.66  eV 

(Fig.  3).  This  value  is  close  to  the  Xe^Cl  energy  formation  in  the  gas 
phase. 


The  experimental  threshold  of  the  Xe2Cl  formation  in  the  Xe  matrix  is 
X  *  370  nm  or  hv  *  3.35  eV,2  so  that  our  theoretical  value  is  overestimated 


by  0.3  eV,  at  least.  It  is  important  to  note  that  the  calculation  of  the 
formation  energy  is  of  lower  accuracy  compared  to  that  of  the  transition 
energy  calculation  [Eq.  (19)],  as  it  involves  more  potentials. 

V .  The  Activated  Xe*Cl  and  Xe^Cl  Complexes  for  Ionic  State  Formation 

As  mentioned  above,  the  Xe^Cl  molecule  is  formed  in  the  Xe  solid  due 

to  the  excitation  of  the  system  by  a  photon.  The  excitation  must  be 

considered  as  the  transition  of  an  electron  to  the  Cl  atom  from  one  of  the 

adjacent  Xe  atoms.  The  ionic  system  which  is  formed  after  the  electron 

transition  (but  before  the  atomic  relaxation)  can  be  considered  as  an 

activated  complex.  According  to  the  Franck-Condon  principle,  the  atoms  of 

the  activated  complex  are  located  at  the  same  points  as  in  the  initial  Cl  + 

Xe  solid  system.  We  shall  consider  two  possible  structures  of  this  complex. 

The  first  is  Xe+Cl  with  the  positive  charge  located  on  one  Xe  atom,  and  the 

second  is  Xe^Cl  with  the  positive  charge  delocalized  among  all  12  Xe  atoms 

around  the  Cl  ion.  Which  of  these  two  structures  is  realized  as  an 

activation  complex  depends  on  the  time  x^  of  the  positive  charge 

delocalization  (or  the  lifetime  of  the  localized  structure).  If  this  time 

is  smaller  than  the  period  of  light  oscillations  x  ,  then  the  Xe^Cl  system 

must  be  considered  as  the  activated  complex,  at  least  for  the  case  where  the 

Cl  atom  is  located  in  the  center  of  the  cage.  For  the  opposite  case  where 

x^j  >>  t^,  the  positive  charge  remains  on  one  Xe  atom  and  the  Xe+Cl  system 

plays  the  role  of  the  activated  complex.  The  lifetime  of  the  localized  Xe+ 

is  on  the  order  of  x^  *  h/2iTEv>  where  e  is  the  energy  of  the  valence 

interaction  between  the  Xe+  and  Xe  atoms.  For  the  underformed  lattice,  the 
+  ° 

Xe  -Xe  distance  is  l  =  4.38  A  and  the  valence  energy  is  approximately  0.3 
eV,^  which  gives  x^  ■  2  x  10  ^  s.  This  time  is  on  the  order  of  visible 
light  oscillations  which  excite  the  Cl  +  Xe  solid  system.  Since  it  is  not 


clear  which  of  the  two  systems  is  the  activated  complex,  we  shall  consider 
both  of  them. 


A.  The  activated  Xe+Cl  Complex 

The  energy  of  the  Xe+Cl  activation  complex  depends  on  the  location  of 
the  Cl  atom  in  the  cage.  The  size  of  the  Cl  atom  is  smaller  than  the 
minimal  Xe-Xe  distance,  so  that  the  Cl  atom  can  be  shifted  in  the  cage 
toward  one  of  the  surrounding  Xe  atoms.  While  these  Xe  atoms  can  also  be 
shifted  toward  the  Cl  atom,  it  is  possible  to  neglect  this  shift,  as  it  was 
shown  to  be  small  (Section  IV. A).  In  order  to  find  the  position  of  the  Cl 
atom  inside  a  cage,  the  potential  of  the  atom  was  calculated  as  a  function 
of  the  shift  6  from  the  center.  It  was  found  that  this  potential  changes 

O 

within  the  limit  of  0.0007  eV  only  for  6  <  0.8  A  and  begins  to  rise 

O 

significantly  after  the  point  6  =  0.9  A.  Similar  results  were  obtained  also 
in  Ref.  7,  where  for  the  temperature  of  12  K  the  thermal  energy  is  kT  * 

0.001  eV,  so  that  the  Cl  atom  can  be  moving  practically  free  around  the  cage 

«  +  . 
center  within  the  boundary  6  <  0.8  A.  We  shall  calculate  below  the  Xe  Cl 

activation  energy  for  both  the  central  Cl  position  (6  *  0)  and  the  shifted 

o 

one  with  6  =  0.8  A. 

The  energy  of  the  Xe+Cl  system  formation  in  solid  phase  with  Cl  at  the 
center  (6=0)  is 


(s) 


vx.*cr(t) 


(8) 


Xe 


*Cl  +  UX.+Cl-(t)  + 


ll0XeCfU) 


12UXeCl(t)  +  V 


(23) 


where  I^e  *  12.13  eV  is  the  Xe  ionization  potential,  A^  *  3.68  eV  is  the  Cl 

(g)  +  . 

electron  affinity,  and  *s  the  Xe  -Cl  potential.  The  polarization 

energy  in  (23)  can  be  calculated  by  using  the  expressions  (5)-(9).  It  is 
found  to  be  equal  to  e  *  -0.49  eV.  The  Xe+Cl  potential  is  equal  to 


^Xe+Cl”'t'  *  '3*37  eV.  Inserting  these  numbers  into  Eq.  (23)  along  with  the 

values  for  \eQl'U)  and  U XeC1<l)  (  see  the  previous  section),  we  obtain 
(s) 

vXe+cfU)  -  4-75  eV  (Fig.  3). 

O 

For  the  case  of  Cl  atom  shifted  by  0.8  A  from  the  cage  center,  all 

interactions  are  the  same  as  for  the  6*0  case  except  for  the  Xe+Cl 

interaction,  which  is  ■  "4.08  eV.  Consequently,  the  activation 

(s) 

energy  is  ^e+Cl"  *  3.94  eV.  Comparing  this  value  with  the  upper  limit 
obatained  for  6*0,  one  finds  that  the  Xe+Cl  activation  energy  in  the  Xe 

solid  matrix  forms  a  band  which  is  in  the  range  3.94  eV  <  <  4.75  eV. 


B.  The  Xe^Cl  Complex. 


Let  us  now  consider  the  Xe^Cl  system  formed  by  the  Cl  ion  located  in 
the  center  of  the  cage  formed  by  12  Xe  atoms  with  the  common  charge  +1.  The 
atomic  Xe  polarization  time  is  much  shoter  than  the  time  of  delocalization, 
and  therefore  the  polarization  of  the  lattice  follows  the  motion  of  the 
positive  charge  among  the  12  atoms  of  the  Xe^Cl  system.  Consequaetly ,  the 
polarization  energy  is  determined  by  the  Cl  ion  and  one  Xe+  ion  and  is 
equal  to  the  polarization  energy  of  the  Xe+Cl  system.  It  follows  that  the 
Xe^Cl  energy  differs  from  that  of  the  Xe+Cl  complex  by  the  valence  energy 
of  delocalization, 


vxe|2cr(t)  ■  V\:f(l)  +  cv(t)  •  (: 

Quantum  chemistry  considerations  show  that  in  the  fee  lattice  the 
delocalization  energy  is 

e  *  /ng  ,  8  <  0  ,  C 

v  o  o 


where  n  is  the  number  of  interacting  neighbor  atoms  and  Ba  is  the  exchange 
integral  between  the  5p  AO's  of  the  two  Xe  atoms  directed  along  the  i-j 


bond.  The  simple  expression  (25)  is  valid  only  for  small  exchange  integrals 


(5^,  which  is  just  the  case  for  Xe  atoms.  Using  the  simple  MO  model  for  the 
Xe^  molecule,  one  obtains 


0  1,2.  2_+v 

'  2(  ru  •  V  • 


(26) 


2  2  + 

where  and  are  the  ground  state  and  the  first  excited  potentials  of 
the  Xe^  molecule,  respectively.  Using  the  ab  initio  potentials  without 
spin-orbit  coupling,  which  are  close  to  the  experimental  potentials,^  we 


can  interpolate  the  numerical  values  of  by  means  of  the  exponential 
function 


Pq(R)  =  -172exp[-1.45R]  ,  R  >  3.0  A  .  (27) 

It  is  reasonable  to  assume  that  in  the  Xe^Cl  system  the  positive  charge 
resides  on  only  the  12  Xe  atoms  surrounding  the  Cl  ion,  so  that  each  of 
these  atoms  is  coupled  with  only  four  other  atoms.  Putting  n  =  4  in  the 
expression  (25),  one  obtains  the  energy  of  delocalization  of  the  Xe^Cl 
system  as 


c  (R)  =  2p  (R)  .  (28) 

V  0 

Using  Vxe+d'  =  4.75  eV  and  8y(i)  *  -0.30  eV,  we  find  that  the  energy  of  the 
formation  of  the  symmetric  Xe^Cl  system  is  V^g+  ^-(O  =  4.15  eV  (Fig.  3). 

The  activation  energy  of  the  Xej^Cl  complex  lies  inside  the  band  of  the 
Xe+Cl  activation  complex  energy.  Consequently,  the  excitation  energy  of 
the  ionic  system  forms  a  broad  band,  both  of  whose  lower  and  upper  limits 
(3.94  eV  and  4.75  eV)  are  determined  by  the  Xe+Cl  excitation  (Fig.  4) 

We  have  considered  the  Xe^Cl  system  as  an  activated  complex  with  the 

O  + 

Xe-Cl  distance  £.  =  4.38  A.  As  the  Xe  atoms  in  Xe^Cl  are  attracted  to  the 
Cl  ion,  they  can  be  shifted  towards  the  Cl  ion.  However,  there  are  two 
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alternative  channels  for  this  shift.  The  first,  which  was  considered  in  the 

previous  section,  is  the  shift  of  two  Xe  atoms  with  the  formation  of  the 

Xe^Cl  molecule.  The  energy  of  the  formation  of  this  molecule  in  the  Xe 

(s) 

solid  was  proven  to  be  equal  to  V  =  3.66  eV.  The  second  possible 

channel  is  the  symmetric  shift  of  all  12  Xe  atoms  and  the  formation  of  the 
quasistable  Xe^Cl  system  with  the  Xe-Cl  distance  smaller  than  the  initial 

O 

l  *  4.38  A.  The  energy  of  such  a  system  is 


(s) 


(8) 


VXe|2Cl'(r)  *  XXe  "  AC1  +  UXe+Cl'(r)  +  ll0XeCl"(r)  '  12UXeCl(i) 


+  UY  (r)  -  U„  MU)  +  e  ( r )  +  e(r)  ,  (29) 

Xe. -M  Xe^n  p  v 


where  Uv  u  is  the  interaction  of  the  12  Xe  atoms  with  the  Xe  matrix  and 

AG  ^  2^* 

one  with  another.  According  to  the  results  of  our  calculations  (Fig.  4), 


the  energy  (29)  has  a  minimum  at  the  point  of  r^  =  4.03  A  with  the  energy  of 

stabilization  (relative  to  the  initial  Vv  +„.-(£)]  of  0.54  eV.  The  energy 

Xg  ^  2'/  ^ 

(29)  at  its  minimum  is  Vv  +  „.-(4.03)  =  3.61  eV.  This  value  is  lower  than 

Xg  22^^ 

the  energy  of  the  Xe^Cl  molecule  (3.66  eV).  However,  the  difference  is 
too  small  to  imply  that  the  Xe^Cl  system  is  more  stable  in  the  solid  phase 
than  the  Xe2Cl  molecule. 

The  energy  of  the  transition  in  the  Xe^Cl  system, 


(s) 


,(s) 


X  +  -  s  v  +  -(r  )  -  Vv®'  (r  )  , 

Xe12Cl  XenCl  Vm'  Xe^Cr  tn'  ’ 


(30) 


is  found  to  be  equal  to  3.89  eV  (X  =  319  nm),  which  is  too  far  from  the 
experimentally  observed  X  =  570  nm.2  At  the  same  time,  the  calculation 
performed  in  the  previous  section  confirms  that  Xe2Cl  is  responsible  for 


the  emission  at  X  -  570  nm.  It  follows  that  Xe^Cl"  rather  than  Xe^Cl"  is 
really  the  most  stable  ionic  system  in  the  Xe  solid. 

One  of  the  assumptions  of  the  present  calculation  is  the  neglect  of  the 
valence  interaction  between  the  Xe  atoms  of  a  molecule  and  the  matrix  atoms. 
Using  the  Xe+-Xe  resonance  integral  (27),  it  is  possible  to  estimate  the 
error  resulting  from  this  neglect.  Let  us  consider  the  Cl  ion  and  three 
neighboring  Xe  atoms  with  a  common  +1  charge.  Two  of  these  Xe  atoms  and  Cl 
ion  belong  to  the  molecule  whereas  the  third  Xe  atom  belongs  to  the  matrix 
and  is  located  farther  from  the  Cl  ion.  Taking  into  account  only  the 
Couloumb  Xe+-Cl  interaction,  we  find  the  valence  interaction  with  the  third 
Xe  atom  to  be  on  the  order  of  0.1-0.15  eV.  It  follows  that  ignoring  the 
nolecule-matrix  valence  interaction,  we  have  made  an  error  which  is  not 
negligibly  small;  however,  it  is  also  not  large  enough  to  significantly 
change  the  results  of  the  calculations. 

VI .  Conclusion 

1.  The  energies  of  ionic  xenon-chlorine  systems  in  a  solid  Xe  matrix 
have  been  calculated  in  order  to  study  their  optical  properties  and 
stability.  The  calculation  takes  into  account  the  electrostatic 
polarization  of  the  matrix  and  the  repulsive  and  van  der  Waals  interactions 
between  all  pairs  of  atoms  involved  (Xe-Xe,  Xe-Cl,  Xe-Cl  ).  The  known 
potentials  of  these  interactions  have  been  interpolated  by  analytical 
expressions. 

2..  The  calculation  of  the  energy  levels  of  the  Xe^Cl  ionic  molecule  in 
the  Xe  matrix  strongly  supports  the  hypothesis  that  these  molecules  are 
responsible  for  the  emission  of  A  =  570  nm  light  observed  experimentally . ^ 
The  calculated  wavelength  of  the  Xe*Cl  optical  transition  is  A  -  540  nm. 
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While  the  shift  of  the  transition  wavelength  in  the  solid  matrix  compared  to 
the  gas-phase  value  (A  =  484  nm)  is  mainly  due  to  the  polarization  of  the  Xe 
matrix  by  the  Xe^Cl  molecule,  the  van  der  Waals  interactions  are  also  of 
importance . 

3.  The  calculation  of  the  energy  of  the  activation  complexes  Xe+Cl  and 
Xe^Cl  reveals  that  the  excitation  energy  needed  for  the  radiative 
formation  of  the  Xe^Cl  ionic  molecule  lies  in  an  interval  between  3.94  eV 
and  4.75  eV. 

4.  The  possibility  of  formation  of  an  ionic  complex  Xe^Cl  in  the 
solid  Xe  has  been  investigated.  The  calculated  energy  of  this  complex  is 
slightly  lower  than  the  energy  of  the  Xe^Cl  ionic  molecule.  However,  the 
transition  energy  of  the  Xe^Cl  complex  is  much  higher  than  the 
experimentally  observed  value,  so  that  most  probably  the  Xe^Cl  complex 
is  actually  unstable  relative  to  the  Xe*Cl  molecule. 
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Lowest-lying  potential  energy  states  of  XeCl  (adapted  from  Ref.  2) 
illustrating  the  photoassociative  production  of  electronically  excited 
XeCl  and  subsequent  fluorescence  of  the  excimer  at  350  no  on  the  C  ■*  A 
band.  The  binding  of  the  ground  state  has  been  exaggerated  for 
clarity. 

Diatomics- in-molecules  potential  surfaces  for  Xe^Cl  leading  to  Xe^  +  Cl 
(adapted  from  Ref.  10). 

Schematic  presentation  of  the  energy  levels  and  transitions  in  the 
Xe^Cl  molecule  (see  text).  The  numbers  stand  for  energies  in  eV. 
Energy  levels  of  the  ionic  systems  Xe^Cl  ,  Xe^Cl  (potential  curve), 
and  Xe+Cl  relative  to  the  neutral  Cl  +  Xe  solid  system. 
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